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Abstract. The study of systems with multiple (not necessarily degenerate) metastable states 
presents subtle difficulties from the mathematical point of view related to the variational 
problem that has to be solved in these cases. We introduce the notion of relaxation height 
in a general energy landscape and we prove sufficient conditions which are valid even in 
presence of multiple metastable states. We show how these results can be used to approach 
the problem of multiple metastable states via the use of the modern theories of metastability. 
We finally apply these general results to the Blume-Capel model for a particular choice of the 
parameters ensuring the existence of two multiple, and not degenerate in energy, metastable 
states. 
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1. Introduction 

In many applicative problems one has to consider a stochastic system evolving in a finite 
not empty state space driven by an energy function. The details of the dynamics depend 
on the system that one has to model, but in most cases the motion is driven by the energy 
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landscape [26]. The typical trajectories are those decreasing the energy; on the other hand 
with small probability the system can perform jumps against the energy drift. 

The smallness of the probability associated to moves against the energy drift is controlled 
by some cost function. This cost is often given by the difference of energy between the two 
states involved in the move. This is the situation that one has to face when studying the 
evolution of Glauber dynamics associated with Statistical Mechanics lattice models [T8] . 
Depending on the model, different cost functions can be introduced. An example, which 
deserved a lot of attention in the recent literature, is that of reversible [21] Probabilistic 
Cellular Automata PEMISEI]. 

The stochastic dynamics is often controlled by a parameter, say the temperature, whose 
value tunes the amount of randomness in the motion by increasing or decreasing the proba- 
bility of the moves against the drift. When the temperature is very low the system freezes 
and tends to stick to the energy drift dynamics ending up to be trapped in the minima of 
the energy. When the temperature is very high, on the other hand, jumps against the drift 
are highly probable and the system moves almost freely in the state space. 

It is then clear that a huge amount of information about the low temperature dynamics 
is obtained when one knows (i) the structure of the absolute minima of the energy (ground 
states) and (ii) the maximal barrier that has to be overcome to reach the set of the ground 
states starting from any possible state, namely, the relaxation height. The severity of prob- 
lem (i) strongly depends on the particular expression of the energy function, but in most 
applications it is not a particularly difficult task. Problem (ii), on the other hand, is always 
a very laborious (and sometimes hard) one since the solution of several variational problems 
is involved. 

The relaxation height is of basic importance in the study of this kind of dynamics. 
For instance, when dealing with the description of the metastable behavior of Statistical 
Mechanics systems, the relaxation height is the quantity controlling the typical time the 
system waits before nucleating the stable state starting from the metastable one [8||22]. 
The problem of computing the relaxation height in specific model is often very hard, in 
particular when there are several states connected to the ground states via cost barriers 
equal to the relaxation height. This situation is found, for instance, in models with more 
than one metastable state \12\ . 

In this paper we shall discuss some generic properties of the relaxation height and, in 
particular, we shall setup a quite general strategy for its computation. Subsequently we shall 
apply the theory to the very interesting case of the Blume-Capel model [H El [HI [17] for a 
choice of the parameter ensuring the existence of two not degenerate in energy metastable 
states. Due to the presence of multiple metastable states this case is particularly interesting 
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and, a priori, very difficult to be treated. 

The Blume-Capel model has been originally introduced in [HHI] to study some particular 
magnetic systems. The model was then generalized to the so called Blume-Emery-Griffiths 
model [5] in order to describe the A transition in He 3 -He 4 mixtures. Those models have 
been widely studied in the literature and many applications have been considered. Their 
particular interest is due to fact that in these three state spin systems on lattice different 
energies are associated with different interfaces. This is not the case in other multi-state 
spin systems, such as the well known Potts model, where a ferromagnetic coupling favors 
the presence of neighboring homologous spins and the same energy cost is assigned to any 
pair of differing neighboring spins. 

This fact is very peculiar of the Blume-Capel model and gives rise to very interesting 
phenomena when the metastable behavior of the system is considered. Indeed it is seen that 
even if one of the three states is not favored from the energetic point of view, it can (and 
indeed does) behave as a bearing between the other two phases if the interfaces between this 
particular state and the other two ones are energetically less expensive (see, for instance, the 
proof of item |2] in Lemma [4.91 where this effect appears in all its importance). 

In [T7] the authors studied metastability in the Blume-Capel model in the so called 
Wentzell-Friedlin regime, that is they considered finite volume, finite positive magnetic 
field, and let the temperature tend to zero. In that paper a wide region of the space of 
parameters of the model was analyzed. More precisely, denoted by A the chemical potential 
(see equation (14. 11) below) and by h the magnetic field, the region h > A > was studied in 
detail. It was proven that in that region the Hamiltonian of the system has a single ground 
state u, namely, the configuration with all the spins equal to plus one, and that the system 
admits a single metastable state, that is the configuration d with all the spins equal to minus 
one. It was also proven that the configuration 0, with all the spins equal to zero, plays an 
important role. For 2A > h > A the transition from the metastable state d to the stable 
one u is direct, although the plus droplet growing inside the minus phase is "protected" by 
a thin frame of zeros. For h > 2A, on the other hand, during the transition the system visits 
the intermediate (not metastable) state 0. This kind of phenomena is not expected to be 
observed in multi-state spin systems such as the Potts model where all the interfaces pay 
the same energy cost. 

In [23] the study of the metastable behavior of the Blume-Capel model was addressed 
in infinite volume in the regime characterized by finite A and h and temperature tending to 
zero. In that paper the authors studied the regions < —A < h and < A < h of the space 
of parameters, see [231 Theorem 1], and proved that, starting from d, in the first one the 
state is observed before the transition to u, while in the second one it is not. 
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Figure 1.1: Schematic description of the energy landscape for the Blume-Capel model ( 14. ip 
at A = and h > 0. The quantity T c represents the energy barrier that must be overcome 
to exit the metastable states. In the picture it is remarked that the two metastable states d 
and are not degenerate in energy. 

This result, looked at from the finite volume point of view, suggests that at negative A 
the state is metastable while at positive A it is not. The first of these two remarks was 
discussed, on heuristics grounds, in [16], while the second was indeed proven in [17]. We can 
then conjecture that in the peculiar case A = and h > both d and are metastable. In 
other words we expect that for this choice of the parameters the Blume-Capel model shows 
the very interesting phenomenon of multiple not degenerate in energy metastable states (see 
figure 11.11) . 

Situations of this kind have been already studied in the case of Probabilistic Cellular 
Automata p2], but there an old approach [25] to metastability had been used. In that 
paper the model dependent study had to be very detailed and, hence, very difficult. In 
this paper we show that by using the more recent approach to metastability, that will be 
shortly reviewed at the beginning of Section H] and the general result on relaxation height in 
Theorem 12.41 below, the problem can be solved relying on few model dependent properties. 

The paper is organized as follows. The general results are discussed in Section [2j Their 
application to the problem of metastability in Statistical Mechanics systems is discussed in 
Section [31 In particular, their application to the study of the metastable behavior of the 
Blume-Capel model is discussed in Section HI Finally we prove the theorems and the lemmas 
in Sections [5] and El 

2. Maximal stability level 

In this section we shall discuss some general results related to energy landscapes. Under very 
general assumptions on the structure of the energy landscape we shall introduce the notion 
of maximal stability level and prove some results yielding a handy recipe for its computation. 
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2.1. Energy landscape 

An energy landscape is a quaternion (X, Q, H, A) where the finite not empty set X, Q C 
X x X, H : X — > IR, and A : Q IR+ are respectively called state space, connectivity 
relation, energy, and cost function, and the relation Q on X is such that for any x, y G X 
there exist an integer n > 2 and xi, . . . , x n G X such that x\ = x, x n = y, and (x;, Xj+i) G Q 
for any % — 1, . . . , n — 1. 

An energy landscape (X, Q, if, A) is called reversible if and only if the connectivity rela- 
tion Q is symmetric and 

H(x) + A(x, y) = A(y, x) + H(y) (2.1) 

for all (x, y) G Q. 
S.S. Definitions 

Consider a reversible energy landscape (X, Q, H, A). Given Y G X such that if (y) = H(y') 
for any y, G F, we shall denote by H(Y) the energy of the states in Y . For any Y G X we 
shall denote by F(Y) the set of the minima of the energy inside Y, that is to say y G F{Y) 
if and only if H(y') > H(y) for any y' G Y. We let X s := F(X) be the set of ground states 
of if, namely, the set of the absolute minima of the energy. 

For any positive integer n, oj G X n such that (tc^Wj+i) G Q for alH = 1, . . . , n — 1 is 
called a pa£/i joining wo to u n ; we also say that n is the length of the path. For any path oj 
of length n, we let 

$ w := max if^) + Afa, 0J i+1 ) (2.2) 

i=l,...,n— 1 

be the height of the patfj. For any y, z G X we denote by f2(y, 2) the set of the paths joining 
y to z and define the communication height between y and z as 

:= min $ w (2.3) 

ti;Gf2(j/,z) 

From (ED, 02211, and Q23J) it follows immediately that 

= $(z,y) (2.4) 

for all y, 2 G X. For any Y, Z C X we let 

<3>(Y, Z):= min <1> = min <fr(y,z) (2.5) 

where we have used the notation Q(Y, Z) for the set of paths joining a state in Y to a state 
in Z. 



Since the energy landscape is reversible, the energy of the state w„ is implicitly taken into account in 
, indeed (|2.ip implies H(uj n ) < A(cj„_i, cj„) + 
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We say that X is fully attracted by X s if and only if X B ) —H{x) = for any x G X\X S 
or (trivial case) X = X s . 

For any x G X we denote by I x the set of states y G X such that < H(x). Note 

that /a; = if x G X s . We then define the stability level of any x E X \ X s 

V x :=$(x,I x )-H(x)>0 (2.6) 

Note that the stability level V x of x is the minimal cost that, starting from x, has to be 
payed in order to reach states at energy lower than H(x). Following [22] we now introduce 
the notion of maximal stability level. 

Definition 2.1 Consider a reversible energy landscape (X,Q,H,A). Assume X \ X s ^ 0, 
we let the maximal stability level be 

T m := sup V x (2.7) 

x<=X\X s 

We also set 

X m :={xeX\X s : V x = T} (2.8) 

Note that, since the state space is finite, the maximal stability level T m is a finite number. 
Note, also, that if X is fully attracted by X s and X\X S ^ 0, then r m = and X m = X\X S . 

2. 3. Results 

The notion of maximal stability level has been introduced by looking at the paths starting 
from any state of X and reaching lower energy states. This point of view is often very useful 
when dealing with metastability problems [22], indeed the maximal stability level controls 
the asymptotic of the time the system spends in the metastable state before nucleating the 
stable one. It is worth noting, on the other hand, that computing the maximal stability level 
of a concrete model is a hard task. Indeed, one has to solve the variational problem $(x, I x ) 
for any x G X. 

It is then clear the interest of results providing sufficient conditions, whose verification 
in the context of specific model is of reasonable difficulty, ensuring that a real number is 
the maximal stability level. We note that this question has already been debated in the 
pertaining literature, see for instance [22l Section 4.2], [H[ Theorem 2.3], [T9], Lemma 1.2 
and hypothesis H2], and [201 Lemmas 1.16, 1.16 and 1.17]. We remark that in all the quoted 
references the authors always stated results in the case \X m \ = 1. 

These types of results can be obtained by looking at the problem by a different point of 
view, that is by looking at all the paths connecting any state not belonging to X s to the set 
of ground states itself. 
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Theorem 2.2 Consider a reversible energy landscape (X,Q,H,A). Assume X\X S ^ 0. If 
the not empty set A C X \ X s and the positive real number a £ IR + are such that 

1. X s ) - H(x) = a for all x £ A; 

2. either X\(AU X s ) = or X s ) - H(x) < a for all x £ X \ (AU X s ); 
then 

T m = a and X m = A 

The above theorem is very useful in the applications, indeed it gives a general strategy 
to approach the problem of computing the maximal stability level in special models. The 
idea is that one has to figure out what is the set of configurations such that starting from 
them the cost to reach the ground states is precisely the maximal stability level. Once it 
has been proven that starting from all the configurations in this set the cost to be payed is 
the same, then one is just left with the proof that starting from any other configuration the 
cost is strictly smaller. And this is not a terrific computation since only an upper bound to 
the height along the paths is needed. 

It is possible to prove a necessary condition in the spirit of the statement in Theorem 12.21 
In other words on the basis of the Definition 12.11 we can say what is the barrier that must 
be overcome to visit the set of ground states X s starting from states in X m . 

Theorem 2.3 Consider a reversible energy landscape (X,Q,H,A). Assume X \ X s ^ 0. 
Then 

1. X s ) - H(x) = T m for all x £ X m ; 

2. X s ) - H(x) < T m for allxeX\ (X m U X s ). 

This theorem ensures that r m is the maximal height that must be overcome, starting 
from any state in the system, to reach the set X s of ground states of the system. For this 
reason the quantity r m will be also called the relaxation height of the system. 

In many situations it is of great help combining the relaxation height and the stability 
level point of view to get sufficient conditions whose verification in the context of specific 
model is of reasonable difficulty We then state the following theorem. 

Theorem 2.4 Consider a reversible energy landscape (X, Q, H, A). Assume X \ X S ^ 0. // 
the not empty set A C X \ X s and the positive real number a £ IR + are such that 

1. X s ) - H(x) = a for all x £ A; 
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2. either X\{AU X s ) = or V x < a for all x E X\(AU X s ); 



then 



r 



in 



a and X, 



in 



.4 



We remark that Theorem 12.21 is a corollary of Theorem I2.4L Indeed, if the set A and the 
real number a satisfy the hypotheses of Theorem [272] then they also satisfy the hypotheses of 
Theorem 12.41 Nevertheless, since in the proof of Theorem 12.41 the hypothesis that the state 
space X is finite play a crucial role while in the proof of Theorem 12.21 it does not, in Section [5] 
the two statements will be proven independently. The idea in the proof of Theorem 12.21 could 
result useful in the study of metastability of Statistical Mechanics systems in infinite volume. 

3. Metastable states of Statistical Mechanics lattice models 

The theory developed above can be fruitfully applied to study metastability in Statistical 
Mechanics lattice models when multiple metastable states are present. A natural setup in 
which this problem can be approached is that of Markov chains or Markov processes. In this 
context powerful theories [Hl|22] have been developed with the aim of finding answers valid 
with maximal generality and to reduce to a minimum the number of model dependent inputs 
necessary to describe the metastable behavior of the system. In this section we shall briefly 
review the pathwise and the potential theoretic point of views and explain, in this context, the 
interest of the results stated in Section [2] in order to deal with the case of multiple metastable 
states. In particular in Section 13.21 we shall prove a recipe to construct metastable sets in 
the potential theoretic approach sense. 

Consider a finite state space X and a function q : X x X — > [0, 1] called connectivity 
matrix such that for all x, y E X 



Moreover assume that for any x, y E X there exist an integer n > 2 and x\, . . . , x n E X such 
that X\ = x, x n = y, and q(xi,Xi+i) > for any i — 1, . . . , n — 1. 

Define Q C X x X by letting (x, y) G Q if and only if q(x, y) > 0. Consider the 
Hamiltonian H : X — > R and the cost function A : Q — > R + such that for all (x, y) G Q 



Moreover assume that the following is satisfied: (i) a state x E X is such that q(x,x) = if 
and only if A(x, y) = for all y G X such that y ^ x and (x, y) G Q; (ii) for all x G X such 
that q(x,x) > the cost function is such that A(x,x) = 0. 




q(x,w) = l and q(x,y) = q(y,x) 



(3.1) 



wex\{x} 



H(x) + A(x,y) = A(y,x) + H(y) 



(3.2) 



metadegenere.tex - 4 marzo 2013 



8 



21:47 



Given 1//3 > 0, called temperature, consider the aperiodic ergodic Markov chain on the 
finite state space X and transition matrix pp : X x X — > [0, 1] given by 



Pp(x,y) := q(x,y)e P A ( x ' y *> Vi,t/ 6 X and x ^ y (3.3) 

and 

p/j(x,x):=l- Pp( x >y) Vx G X (3.4) 

From the second equality in (13. ip and from (13.21) it follows that the chain is reversible with 
respect to the Gibbs measure 

4- e-W*) (3.5) 

where is the partition function, that is to say iJLp(x)pp(x,y) = pp(y,x)[ip(y) for any 
x,y G X. 

Noted that the quaternion (X, Q, H, A) is a reversible energy landscape, see Section |2~T| 
we use the theoretical setup introduced in Section [2] and, following [22] , we call the set X m 
introduced in Definition 12.11 set of metastable states. In the following two subsections we 
shall recall why such a set deserves its name. 

3.1. Pathwise approach 

The reason way X m deserves its name has been widely explained in [22] in the framework 
of the so called pathwise approach to metastability. In that paper the important properties 
of the states in X m have been fully studied in the case of the Metropolis dynamics, namely, 
when the cost function is chosen as A(x,y) := [H(y) — H(x)] + , where, for any real a, [a] + is 
equal to a if a > and to otherwise. Note that this choice for the function A is coherent 
with the theoretical setup introduced above since ( 13 .2\\ is satisfied. 

A partial generalization of the results in [22] to the present case, that is to say when 
the function A is as general as explained in Section [21 has been done in [13]. The authors 
had to consider this more general situation in [14] since such a setup arises naturally when 
Probabilistic Cellular Automata are studied. 

Here, in order to justify the name of set of metastable states given to X m , we just recall 
the following result proven in [T3]: for any e > and x G X m 

lim P x (e /3(rm - e) < t Xs < e^ Tm+e) ) = 1 (3.6) 



where F x is the probability measure on the space of trajectories of the Markov chain started 
at x and Tx B is the first hitting time to X s for the chain started at x. In the case of 
the Metropolis dynamics, this result has been proven in [22j Theorem 4.1]; in that paper a 
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much more detailed description of the metastable behavior of Metropolis dynamics, including 
results in distribution and in law [22j Theorems 4.9 and 4.15], has been given. 

One of the remarkable features of the pathwise approach is that a constructive definition of 
metastable states is given. When studying a particular model, in order to find the metastable 
states, one just has to find the set X m , that is one has to compute the maximum stability 
level T m and identify all the states whose stability level is equal to T m . 

Unfortunately this is in practice a laborious task. Indeed the solution of several varia- 
tional problems is required. It is then very useful the strategy suggested by the Theorems 12.21 
and 12.41 for the computation of r m . This approach has been already used in some situations 
but always in cases in which there existed a single metastable state, see for instance j22j Sec- 
tion 4.2] and [H Theorem 2.3]. 

3.2. Potential theoretic approach 

Different interesting approaches to metastability have been developed in the recent literature, 
see for instance [2ll8HT0ll20] . The notion of metastable states is given in different ways; but 
these different notions are indeed strictly related. For instance in [8] it is introduced the 
"metastable set" which, by using the language introduced above, can be reinterpreted in 
some cases as X m U X s . We discuss this issue in this subsection in detail. 

Following [TJIH] we apply the theory developed in the Section [2] to the model defined at 
the beginning of Section [3j The potential theoretic approach to metastability, introduced 
in [7J[8], has been lately developed for Kawasaki dynamics in [9j[20] and for probabilistic 
cellular automata in |24j . 

In this approach the main definition of interest is that of metastable set [8] . For a generic 
Markov chain the Dirichlet form is defined as the functional 

£fi[h] = \ £ VLp(x)p p {x,y)[h(x)-h(y)]* = ± £ ^e-^ + ^[h(x) - h(y)} 2 (3.7) 

x,yeX x,yeX @ 

where h : X — > R is a generic function and we have used (13. 3p and the definition (13. 5p of 
Gibbs measure. 

Given two non-empty disjoint sets A, B C X the capacity of the pair A and B is defined 

as 

CAPb(A,B):= min £ p (h) (3.8) 

Note that the capacity is a symmetric function of the sets A and B. It can be proven that 
the right hand side of (13. 8 p has a unique minimizer called equilibrium potential of the pair A 
and B and given by 

h* A B {x) = F x (t a < t b ) 
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for any where ta and tb are, respectively, the first hitting time to A and B for the 

chain started at x. 



Definition 3.5 A set M C X is said to be p.t.a.-metastable if 

lim max^ M /x^(x)[CAPff(a;, Af)]" 1 = Q 
/3™ min^M/i^^pAP^^M \ {a;})]" 1 

The prefix p.t.a. stands for potential theoretic approach. We used this expression in order 
to avoid confusion with the set of metastable states X m introduced in Definition 12.11 

In [5] and in the related papers the properties of the p.t.a.-metastable sets have been 
widely described. Here we just mention the following one to explain why these states are 
called metastable; we refer the interested reader to [HE]. Let M be a p.t.a.-metastable set. 
Let x e M and J C M \ {x} be such that for all y G M \ ( J U {x}) either 

OE-^J- < 1 or py ' < 1 

^(^) CAP^^J) ^ 

where by we mean that the ratios on the left-hand sides are uniformly bounded from 
above by a function of (3 tending to zero in the limit /3 — > oo. Then 

Ei[tj1 = capm (1 + o(1)) 

and 

F x (tj > m x rj) = [1 + o(l)]e-*[ 1+0 W] t > 

for /3 large enough, where is the average for the chain started at x and, for any x e M, 
the valley A(x) C X is defined as 

A(a;) := {yeX: ¥ y (r x = t m ) = sup F y (r z = r M )} 

When the potential theoretic approach is applied to study metastability in specific models 
the main difficulties are the identification of the p.t.a.-metastable sets (e.g., the set made 
by the two configurations with all the spins respectively equal to minus one and to plus one 
in the standard Ising model) and that of giving sharp estimates to the capacities among the 
elements of this set. In this subsection we will show how the first of the two points above 
can be approached by using the notion of metastable states introduced in Definition 12 . 1 L 

It is important to note that the Definition 12.11 of set of metastable states and the Defi- 
nition |33] of p.t.a.-metastable sets are different in spirit. The set of metastable states X m 
is defined univocally, that is to say that all the states that are thought to deserve the name 
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of metastable states are collected in the same set X m . In other words Definition 12.11 is con- 
structive and the set of metastable states X m is unique, even though its cardinality can be 
larger than one. On the other hand the p.t.a.-metastable sets are defined as those sets of 
states satisfying condition (13.91) . so that several p.t.a.-metastable sets can exist. This sub- 
tle, but important, difference is often hidden when systems with a unique metastable state 
(\X m \ = 1) are considered, but it emerges in all its importance when multiple metastable 
states (\X m \ > 2) are present. 

This aspect of the potential theoretic approach is, indeed, one of its distinguishing points. 
It reveals the high degree of tunability of the theory. The states in X m are metastable in the 
sense that in order to decrease the energy, starting from them, the highest energy barrier 
(the relaxation height) has to be overcome; this is for sure a definition of metastable states 
absolutely close to the empirical physical meaning of the word metastable. In Theorem 13.61 
we will show how to construct p.t.a.-metastable sets by means of states in X m and X s ; in 
some sense we will construct p.t.a-metastable sets by wisely collecting "maximally stable" 
metastable states. 

One could also enlarge its point of view by considering all the states such that starting 
from them the energy barrier that has to be overcome in order to lower the energy is larger 
than or equal to V, for some V G M such that < V < T m . This set would obviously 
contain X m , but, for V small enough, also some states not in X m will belong to such a set. 
These states are not "maximally stable" metastable states, but, to some extent, they can be 
considered metastable since the energy barrier V has to be payed in order to decrease the 
energy. The potential theoretical approach is well suited to study also these states, indeed 
one just have to look for the p.t.a.-metastable states inside the new set just defined above. 
In the framework of the pathwise approach something in this sense has also been proven 
in [22], see, for instance, the definition [22j equation (2.13)] of metastable set at a prescribed 
level and the related result [22} Theorem 3.1]. 

In this subsection we try to clarify the connection between the "maximally stable" 
metastable sets in the two specified senses and, in particular, we will show how to con- 
struct the p.t.a.-metastable sets once the set of metastable states X m is known. This seems 
to be a very smart recipe since X m can be identified on the basis of the theory developed in 
Section [2J for instance by means of the sufficient conditions given in Theorem 12.21 or those 
in Theorem 12.41 

Define on X the following relation: x, y G X are x ~ y if and only if 3>(x, y) — H(x) < T m 
and x) — H(y) < T m . It is easy to prove that ~ is an equivalence relation. Assume, now, 
that X \ X s 7^ 0, so that X m is not empty. We denote by Xm , . . . ,lm m ' the equivalence 
classes in which X m is partitioned with respect to the relation ~ (see figure |3~2]) . Analogously, 
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m 



Figure 3.2: Schematic representation of the energy landscape aimed to illustrate the defini- 
tion of the equivalence relation ~ and the partition in equivalence classes of the sets X m of 
metastable states and X s of ground states. 

(1) (k ) 

we denote by Xs , . . . ,Xs the equivalence classes in which the relation ~ partitions X s 
(see figure [372]) . 

Theorem 3.6 Assume that X\X S ^ and X\(X s UX m ) ^ 0. Chose arbitrarily x s ^ G X s 
for any i — 1, . . . , k s and x m ^ G Xm for any i — 1, . . . , k m . The set 



is p. t. a. -metastable. 

In words we can say that any set constructed by considering one and only one element 
of each equivalence class Xm and X s is a p.t.a.-metastable set. 

4. Metastable behavior of the Blume-Capel model 

In this section we apply the theory described above to the specific case of the Blume-Capel 
model with zero chemical potential. 

As it has already been pointed out in the pertinent literature, in the study of the 
metastable behavior of a Statistical Mechanics model less details on the structure of the 
energy landscape of the model yield less instructive results on the metastable behavior of 
the system. This remark is quite obvious, but nevertheless interesting, since the study of 
the metastable behavior of a system is typically a very difficult task. So it is useful to un- 
derstand fully to which extent the description of metastability can be pushed forward once 
some model dependent results are known. 

Our results will be then stated according to the following main idea: we shall state 
model dependent lemmas containing informations on the energy landscape of the Blume- 
Capel model followed closely by theorems stating properties of the metastable states of the 
system that can be deduced by using the related lemma. 
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4-1. The Blume-Capel model 

Consider a finite squared torus (periodic boundary condition) A C Z 2 endowed with the 
Euclidean distance d : A x A — > M. + on the torus. As usual we shall misuse the notation and 
let 

d(J, J) : = min d(z,j) 
ieijeJ 

for any J, J C A. We say that two sites i, j G A are nearest neighbor if and only if d(i,j) = 1. 

Let { — 1, 0, +1} be the single spin state space and X := { — 1, 0, +1} A be the configuration 
space. The Hamiltonian of the model is 

G(a) = (*(<) - <3)f - A Yl <^ - h E *(«') 

<i,j> iGA iGA 

for any a £ X, where the first sum is extended to the pairs of nearest neighbors, A G K. 
is called chemical potential, and h G K is called magnetic field. The metastable states for 
/i > A > have been widely studied in [T6|[T7]. 

In this paper we shall attack the case A = which is particularly relevant because two 
metastable states will be proven to exist; in other words from now on we shall consider the 
zero chemical potential Blume-Capel model defined by the Hamiltonian 

H{o)= £(a(i)-a(j)) 2 -/i5>(i) (4.2) 

<i,j> «GA 

Given a G X, H(a) will be also called energy of the configuration a. 

The stochastic version of the model is the Markov chain o~ t , with t = 0, 1, . . . the discrete 
time variable, with transition probabilities (13.31) with A(er, 77) := [H{rj) — H(o~)] + and the 
connectivity matrix defined as 

if <t, 77 differ at more than one site 

—7—7 if a, 77 differ at one single site 
2|A| 

for any a, 77 G X and a ^ 77, and q(o~, a) = 1 for any er G Af. Note that, since the dynamics 
is Metropolis, all the results in [22] hold for this model. 

It is worth nothing that from the definition of the transition matrix q given above, it 
follows that for any a, 77 G X there exist a path connecting a to 77, namely, |fi(<7, 77)] > 1. 

4-2. Metastable states of the Blume-Capel model 

We shall discuss the metastable behavior of the zero chemical potential Blume-Capel model 
for the following choice of the parameters. For any positive real a we let \_a\ be the largest 
integer smaller than or equal to a. 
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4-1 



Figure 4.3: Schematic representation of the configuration V c . The configuration Q c has the 
same geometry with the zeros replaced by the pluses and the minuses replaced by the zeros. 
Note that the protuberance can be either on the left or on the right vertical edge (the longest 
ones in the picture) and there it can be shifted freely. 

Condition 4.7 The magnetic field h and the torus A are such that < h < 1, [2/h\ is not 
integer, and |A| > 49//i 4 finite. 

For h > it is immediate to remark that the set of ground states of the energy is 
X s = {u}, where u £ X is the configuration such that u(i) = +1 for all i G A. Indeed the 
exchange interaction, i.e., the positive defined first term of the Hamiltonian in (14. 21) . gives 
its minimal contribution which is equal to zero. And so does the magnetic field part. 

Other two very relevant configurations are d and 0, that is the configuration in which 
all the spin are minus one and the one in which all the spins are zero. Indeed in these 
configurations the exchange part of the energy is minimal, although the magnetic part is 
not. We also note that 

H(u) = -\A\h < H(0) = < H(d) = +\A\h 

On physical grounds, see (for instance) the zero temperature phase diagram in [I7J Fig. 1] , 
it is quite reasonable to guess that d and are possible metastable states of the system. 
This is indeed true as it will be stated in the following Theorem 14.101 To get this result we 
shall use the strategy outlined in Section |2j in particular we will use the Theorem 12.41 with 
guess set A = {d,0}; this will yield the sufficient model dependent statements. 

We now define the critical length of the model as 

4 := [\\ + 1 (4.3) 
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and remark that, by Condition I4.7[ it follows that 

- < 4 < - + 1 and i c > 3 (4.4) 
h h v ' 

We denote by V c the set of configurations in which all the spins are minus excepted those, 
which are zeros, in a rectangle of sides long £ c and £ c — 1 and in a site adjacent to one of the 
longest sides of the rectangle (see figure l4~3l) . We denote by Q c the set of configurations in 
which all the spins are zeros excepted those, which are pluses, in a rectangle of sides long £ c 
and £ c — 1 and in a site adjacent to one of the longest sides of the rectangle (see the caption 
of figure H]3]). From (|4~21) it follows that 

H(V C ) - H(d) = H(Q c ) - H(0) = A£ c - h[£ c (£ c - 1) + 1] (4.5) 

We then set 

r c := H(V C ) - H(d) = H(Q C ) - H(0) (4.6) 

A simple direct computation shows that for h small one has r c ~ 4/h. 

We are now ready to state the model dependent lemmas on which it will be based the 
description of the metastable behavior of the zero chemical potential Blume-Capel model 
(14. 2 p (see Section [6] for the proof). 

Lemma 4.8 Consider the zero chemical potential Blume-Capel model jjj4.2\ ) and assume that 
Condition \4-7\ is satisfied. We have that 

1. for any configuration rj G V c there exists a path u G f2(d, u) such that $^ — H(d) = T c , 
Ui = f] for some i G {l,...,n} (where n is the length of the path), and H(u>j) < 
H(d) + T c for all j G {1, . . . , n} and j ^ i; 

2. for any configuration rj G Q c there exists a path u G f2(0, u) such that $ w — H(0) = T c , 
Ui = rj for some i G {1, . . . , n} (where n is the legth of the path), and H(ujj) < H(0)+T c 
for all j G {1, . . . , n} and j ^ i. 

Lemma 4.9 Consider the zero chemical potential Blume-Capel model ( |^.i?| ) and assume that 
Condition \4-7\ is satisfied. We have that 

1. $(d, u) - H(d) = $(0, u) - H(0) = T c ; 

2. V a < T c for any o G X \ {d, 0, u}. 

Note that the above lemma ensures the validity of the hypotheses of Theorem 12.41 with 
A = {d, 0}. Then, by using the general results discussed in the above section, one can 
identify the set of metastable states and the relaxation height (see Section |6] for the proofs). 
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Theorem 4.10 Consider the zero chemical potential Blume-Capel model UTW and assume 
that Condition \4-7\ is satisfied. We have that 



X m = {d, 0} and T m = T c 

Moreover we can prove a theorem giving the asymptotic of the exit time for the system 
started at the metastable states. 

Theorem 4.11 Consider the zero chemical potential Blume-Capel model (JJty and assume 
that Condition \4-7\ is satisfied. The following holds: 

1. considered the chain started at d, for any e > 



lim P d (e /3(rc " £) < r u < e ^ Tc+e) ) = 1 

and 

lim -E d [r u ] = T c 

/3->oo p 

where r u denotes the first hitting to u of the chain started at d and Pd and Ed denote, 
respectively, the probability and the average on the space of trajectories started at d; 

2. Considered the chain started at 0, for any e > 

lim P (e /3(r ^ e) < r u < e ^ Tc+£) ) = 1 

and 

lim -rEo[r u ] = T c 

/3->oo p 

where, here, r u denotes the first hitting to u of the chain started at and Po and Eo 
denote, respectively, the probability and the average on the space of trajectories started 
at 0. 



4-3. Escape mechanism 

One important result in metastability, besides proving the asymptotic on the exit time, is 
that of identifying the escape mechanism. The typical result is that, in order to perform the 
transition to the stable state, the system has to nucleate a critical droplet of the stable phase 
inside the metastable one. To prove this result supplementary model dependent properties 
are needed. 

Lemma 4.12 Consider the zero chemical potential Blume-Capel model lji4-2\ ) and assume 
that Condition \4-7\ is satisfied. We have that 



metadegenere.tex - 4 marzo 2013 



17 



21:47 



1. any path u G Q(d, u) such that ^ bJ — H{d) = T c visits V c , that is there exists an integer 
i such that u t G V c ; 

2. any path u G f2(0, u) such that $ w — iJ(O) = T c visits Q c , that is there exists an integer 
i such that Ui G Q c . 



Theorem 4.13 Consider the zero chemical potential Blume-Capel model (4JZ) and assume 
that Condition ^. 7| is satisfied. The following holds: 

1. considered the chain started at d, we have that 

lim F d (r Vc < T U ) = 1 

/3->oo 

where r-p c and r u denote respectively the first hitting time to V c and u for the chain 
started at d; 

2. considered the chain started at 0, we have that 

lim P (r Sc < T U ) = 1 

p— too 

where, here, tq c and r u denote respectively the first hitting time to Q c and u for the 
chain started at 0. 



4-4- Remark on the proof of the nucleation property 

When studying the metastable behavior of a system, the detail of the results that one gets 
depends on the amount of model dependent properties that one is able to prove. In order to 
prove the time asymptotic in Theorem 14 . 1 1 1 only Lemma l4~9l is needed. As already remarked, 
to get the nucleation property in Theorem 14.131 it is needed the Lemma 14.121 in which it is 
identified the state that must necessarily be visited by a path joining the metastable state 
to the stable one with maximal height equal to the maximal stability level plus the energy 
of the starting metastable state. In particular it is proven that such a state is in V c for the 
paths starting at d and in Q c for those starting at 0. 

Following the nomenclature in [22], to which we refer the interested reader for a more 
detailed discussion, we say that, given a, t] G X, a set y C X is a gate for the pair of 
configurations a and t] if and only if all the configurations in 3^ have energy equal to the 
communication height between a and t] and any path joining a to 77 with maximal height 
equal to the communication height between a and rj has to pass necessarily through y itself. 
More precisely, a set y C X is a gate for the pair of configurations a and r] if and only if 
H(() = $(cr, 77) for any ( G y and for any path oj G fl(a,rj) such that ^ = 77) there 
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must be ( £ y and % integer such that Ui = (. Thus, Lemma 14.121 implies that the set 
V c C X is a gate for the pair of configurations d and u, whereas the set Q c C X is a gate 
for the pair of configurations and u. 

Less details on the structure of the energy landscape of the model would have yielded 
less instructive results. Indeed, in the proof of Lemma 14. 9[ in order to get the statement we 
shall define the set X d of all the configurations having precisely £ c (£ c — 1) + 1 minus spins. As 
a byproduct of the proof of Lemma 14.91 we shall also get (for free) that any path connecting 
the metastable state d to the stable state u with maximal height equal to the maximal 
stability level plus the energy of the starting metastable state must necessarily visit the set 
F(X d ) of configurations in X d with minimal energy. So that, without any supplementary 
effort with respect to the proof of Lemma [4791 we have that F(X d ) is also a gate for the pair 
of configurations d and u. Then, by the theory developed in [22], we get 

lim F d (r F (x d ) < r u ) = 1 

where the chain is started at d and Tp(x A ) and r u denote respectively the first hitting time 
to F(X d ) and u. 

This result can be looked at as a nucleation statement, but less precise if compared to 
that in the first part of Theorem 14.131 This is somehow natural since it has been obtained 
with less information on the structure of the energy landscape of the system. In the language 
of [22] both V c and F(X d ) are gates for the pair of configurations d and u, but V c is a so 
called minimal gate. 

A similar discussion can be repeated for the pair of configuration and u with X d 
replaced by Xq, namely, the set of configurations in which all the spins are zeros excepted 
for £ c (£ c — 1) + 1 which are pluses. 

4-5. Remarks on gates and minimal gates 

In [221 page 604] a gate y C X for the pair of configurations a and r\ is said to be minimal 
if and only if for any proper subset y of y there exists a path joining a to r\ with maximal 
height equal to the communication height between a and rj which does not pass through y . 
In other words a gate is minimal if and only if all its proper subsets are not gates. In view 
of this we have immediately that the gate X d for the pair d and u and the gate Xq for the 
pair and u are not minimal. On the other hand from Lemma 14.81 we have that V c and Q c 
are minimal gates for the specified pair of configurations. 

In general minimal gates are not unique. In [221 Theorem 5.1] it is proven that the union 
Q (a, rj) of all the minimal gates for the pair of configurations a and rj is the set of what 
the authors call essential saddles (the interested reader is referred to [221 page 604] for the 
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definition). In our model we can prove that V c is the unique minimal gate for the pair d and 
u and Q c is the unique minimal gate for the pair and u, so that we can state the following 
theorem. 

Theorem 4.14 Consider the zero chemical potential Blume-Capel model lj4-^\ ) an d assume 
that Condition \4- 7| is satisfied. Then 

G(d,u) = V c and g(0,u) = Q c 

We finally remark that by [22| Theorem 5.1] it follows that in our model configurations 
in V c (resp. Q c ) are the sole essential saddles for the pair of configurations d and u (resp. 
and u). 



5. Proof of the general results 

In this section we prove the theorems and lemmas concerning the general results stated in 
Section [2] and Section [3] 

Proof of Theorem 12.21 First of all we note that, if X \ (A U X s ) ^ one has that 

V x < a ViGl\(AUl s ) (5.1) 

Indeed, by item [2] in the hypotheses it follows that there exists a path u £ Q(x,X s ) such 
that <3>^ — H(x) < a. Since H(x) > H(X S ), we have that the same path belongs also to 
f2(x, I x ) so that V x < ^ — H(x), which proves (15. ip . 
We want to prove, now, that 

V x = a Vx £ A (5.2) 
The same argument developed above and item [1] in the hypothesis prove that 

V x < a Vx £ A (5.3) 

We are now left with the proof of the lower bound 

V x > a Vx £ A (5.4) 

Pick x £ A and, by absurdity, assume that V x < a. Then, there exists X\ £ X and X\ ^ x 
such that H(xi) < H(x) and a path uji £ Q(x,Xi) such that 

3V - H(x) < a (5.5) 

If it were X\ £ X s we would immediately get a contradiction, since we would have u>\ £ 
Q(x,X s ) and hence $(x, X s ) — H(x) < a in contrast with the item [1] of the hypothesis. 
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Assume, finally, that x\ G X \ X s . By items [T] or [2] in the hypothesis we have that there 
exists a path u 2 G Q(xi,X s ) such that 

^ - H{ Xl ) < a (5.6) 

Now, note that by gluing the path u 2 to the last configuration of u)\ we get the path cjiCJ 2 G 
Q(x,X s ). Moreover, 

$ WiaJ2 = max{$ 1Jl ,$ 1J2 } 

and hence 

" #(z) = max{$ Wl - #(2), $ W2 - if(x)} (5.7) 

Note, now, that 

$ W2 - if (x) = - if (xx) + if (xi) - if (x) 
Since if (xi) < if (x), from (15.61) it follows that 

$„ 2 - if (x) < a (5.8) 

In conclusion, by (15. 5p . (15. TJ) , and (15. 8p we have that 

- if (x) < a 

which, recalling that G Q(x,X s ) contradicts item [1] in the hypothesis. By contradiction 
we have that V x > a and hence prove (15. 4p . 

By (JEl and it eventually follows flO]) . Recalling flU} in Definition O, from (1BT21 
and (15.11) we get the desired equality T m = a. 

We finally remark that, since r m = a, from Definition 12.11 it follows immediately that 
X m = A. Indeed, since T m = a and X m is the collection of states in X \ X s such that 
V x = r m , from items [1] and [2] in the hypothesis it follows that X m = A. □ 

Proof of Theorem 12.31 By Definition 12.11 it follows immediately that 

$(x, 4) - H{x) = T m Vx G X m (5.9) 

and 

$(x,4)-Jf(x) <T m Vx6l\(I m Ul s ) (5.10) 
Proof of item [T] first of all we note that from (15.91) it follows immediately that 

$(x,X s )-if(x) >T m VxGX m (5.11) 

Indeed, recalling (12. 5p . we have that X, C f x implies $(x, J x ) — if (x) < $(x,X s ) — if(x). 
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To prove the upper bound, we pick x G X m and note that, by (15. 9p . (I5.10p . and 
the finiteness of the state space X, we have that there exist a sequence of n — 1 states 
X\, . . . , x n -\ G X\X 8 , a ground state x n G X s , and n paths u>i G f2(xj_i, Xj) with z = 1, . . . , n 
such that 

H(x ) > H( Xl ) > H(x 2 ) >■■ H{x n ^) > H(x n ) 

and 

$ Wi - if (a;,.!) < r m 

for all z = 1, . . . , n, where we have let xo = x. Note that the bound is not strict because 
some of the states x , X\, . . . , x n _i could belong to X m (at least x = x does belong to X m ). 
We then remark that the path U1U2 ■ ■ ■ ou n obtained by gluing the above paths belongs to 
Q(x, X s ) and 

Qawa-un ~ H{x) = max [$ Wi - H(x)] (5.12) 

i=l,...,n 

For any z = 2, . . . , n we have that 

- if (x) < r m and ^ - if (s) = $ Wi - ff (xi-i) + if (x^x) - if (X) < T m (5.13) 

since if (xj_i) < X(x) for all % — 1, . . . , n. In conclusion we have that 

$ _ H(x) < r 

which, recalling that oj\uj<i • • • 0J n G fi(x, X s ) proves the upper bound 

$(x, X s ) - H(x) < T m Vx G X m (5.14) 

Equations (15. lip and (I5.14p finally yields item [1] of the theorem. 

Proof of item [2j we follow the same strategy as that used in the proof of the upper bound 
(I5.14p above. We just notice that, since x G X\ (X m UX s ), in equation (15.131) the first upper 
bound is strict, as it follows from (I5.10p . □ 

Proof of Theorem 12.41 If X \ (A U X s ) 7^ 0, by following the same strategy described in 
the proof of Theorem 12.21 one easily gets 

V x < a Vx G A (5.15) 

Since one has to prove that V x = a for any x G A, a lower bound is needed. Then, given 
a generic x G A, we assume by absurdity that V x < a. Thus, there exists x\ G X such that 
if(xi) < if(x) and a path lo\ G fi(x,Xi) suc h that 

$ W1 - if (x) < a (5.16) 
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If it were X\ G X s we would immediately get a contradiction, since we would have 
oji G fi(x,X s ) and hence <&(x,X s ) — H(x) < a in contrast with item [1] in the hypotheses of 
this theorem. 

If, on the other hand, x\ G X \ X s , by exploiting the fact that X s is a finite set, we can 
proceed as above and find a sequence of n — 2 states x 2 , ■ ■ ■ , £ n -i 6l \ X,, a ground state 
x n G X,, and n — 1 paths G Xj) with i — 2, . . . , n such that 

H( Xl ) > H(x 2 ) > ■■■ F(x n _i) > #(x„) 

and 

$ Wi - H{xi- X ) < a 

for all % = 2, . . . , n. Note that the bound is not strict because some of the configurations 
xi, . . . , x n _i could belong to A. We then remark that the path Uiu 2 ■ ■ • co n obtained by gluing 
the above paths belongs to Q(x,X s ) and 

Qwvn - H{x) = max - H{x)\ (5.17) 

«=l,...,n 

Let Xq = x, for any i = 1, . . . , n we have that 

$ Wi - H(x) = $ Wi - H(x i - 1 ) + H{xi- X ) - H(x) < a 
since if < X(x) for all % — 1, . . . , n. In conclusion we have that 

$wiw 2 -w„ - -f^(x) < a 

which, recalling that W!W 2 • • • w n G £l(x, X s ) contradicts item [1] in the hypotheses of this 
theorem. 

We have finally proven that V x = a for any x G A. This, together with item [2] in the 
hypotheses of the theorem, implies that r m = a. 

We finally remark that, since r m = a, from Definition 12.11 it follows immediately that 
X m = A. Indeed, since T m = a and X m is the collection of states in X \ X s such that 
V x = T m , from items [T] and [2] in the hypothesis it follows that X m = A. □ 

We consider now the model in Subsection 13.21 and prove the Theorem 13.61 Before prov- 
ing the theorem we recall the following elementary estimate given in [9] in the context of 
Kawasaki dynamics. 

Lemma 5.15 For every non-empty disjoint sets A,BcX, there exist constants < C\ < 
C*2 < oo ( depending on A and B ) such that for all /3 > 

d < e mA ' B) Zp CAP^A B) < C 2 (5.18) 
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Proof of Lemma 15.151 the lemma can be achieved via the same argument developed in 
the proof of [9j Lemma 3.1.1]. □ 



Proof of Theorem 13.61 Let M := {x Sj i, . . . , x St k s , x m> i, . . . , x m ,fc m }- Strategy of the proof: 
we bound the numerator in the left-hand side of (13.91) from above and the denominator from 
below. These bounds will be sufficient to prove that the ratio tends to zero for /3 — > oo. 

We first bound the numerator in the left-hand side of (13.91) from above. By using the 
lower bound in (15 . 1 8[) we have that 

^(x)[CAP p(x,M)r x < to(x)±.e?*<* u >Z, l = -L g/WMO-HM] (5 19) 

for any x M. 

Let x G X\(X m UX s ). Since Mnl s W ^ for any % = 1, . . . , fc B , item [2] of Theorem O and 
the definition of the equivalence relation ~ given in Section [321 imply that <$>(x, M) — H(x) < 

r m . 

Let, now, x G (X s UX m )\M (note that this set can be possibly empty). Since 
for any i = 1, . . . , fc s , and M H Xm ^ for any i = 1, . . . , fc m , from the definition of the 
relation ~ given in Section I3T2"1 we have that $(x, M) — H(x) < T m . 

These remarks and (I5.19P yield that there exist two real numbers C3 < 00 and 5 > 
such that 

fip(x)[CAP p (x, M)]- 1 < C 3 eP (Tm ~ 5) (5.20) 

for any x ^ M. 

We, now, bound the denominator in the left-hand side of (13 .9p from below. By using the 
upper bound in (I5.18P we have that for any x G M 

^(xjlCAP^x, M \ {x})]^ 1 > ^(^^^(^^^Z, = i ^W- M \W)^W] (5 . 21) 

First assume x G M(lX m . From the definition of the equivalence relation ~ and the fact 
that \M n Xm\ = 1 for any i — 1, . . . , k m , we have that, provided \M D X m \ > 2 , 

<5>(x,(M\{x})nX m )-H(x)>T m 

On the other hand, since M D Xs 7^ for any % = l,...,k m , from item [1] in Theorem 12.31 
and from the definition of the equivalence relation ~ we have that 

$(x,(M\{x})r\x B )-H(x) >r m 

Recalling (I2.5p . from the two inequalities above, we have that 

M \ {x}) - H{x) = min $(x,y) - H(x) > T m 

yeM\{x} 
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for any x G M PI X m . 

Suppose, now, x G M n X s . By proceeding as above, recalling (12. 4p . and using that 
H(X S ) < H(y) for any y G X m , we prove that 

$(x,M\ {x}) - if (a;) > r m 

for any x G M D X s . 

We have eventually proven that $(x, M \ {a;}) — if (re) > r m for any x G M. This remark 
and (I5.2ip yields the lower bound 

MxHCAP^MUx})]- 1 > (5.22) 

for any x G M. Hence, from (15.201) and (I5.22p we finally get 

max^ M /f^(x)[CAP^(a;,M)]~ 1 C3 _^ 
min x6M ^(x)[CAP^(x, M \ {x})]- 1 ~ C 2 6 

that completes the proof of the theorem. □ 



6. Proof of results concerning the Blume— Capel model 

In this subsection we prove the Lemmas 14.81 and 14.91 and the Theorems 14.101 and 14.111 Before 
starting the proof of Lemmas 14.81 and 14.91 we state few technical results on polyominoes based 
on the paper [1] to which we refer the interested reader for more details. 

6.1. Some results on polyominoes 

Polyominoes have been widely studied both by physicists, to model crystal growth, and by 
combinatorialists. In combinatorics the main problem has been that of counting the number 
of polyominoes given their area or their perimeter [6]. 

In the framework of metastability polyominoes play an important role and the typical 
relevant problem is that of finding the polyominoes with minimal perimeter and given area. 
This problem has been widely studied by Alonso and Cerf [1] both in dimension two and 
three. In the following lemmas and in the corollary, we shall summarize some of the properties 
proven in that paper that we shall use in the sequel. 

In order to state the lemma we need some preliminary definitions. Consider the lattice 
Z 2 embedded in R 2 . Two sites of Z 2 are said to be nearest neighbor if and only if their 
mutual Euclidean distance is equal to one. A unit square is a square of area one, whose 
center belongs to Z 2 and whose vertices belong to its dual Z 2 + (1/2, 1/2). A polyomino is 
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Figure 6.4: From the left to the right: two not connected convex polyominoes, a connected 
not convex polyomino, a connected convex polyomino which is not a convex subset of M 2 , 
a connected convex polyomino which is a convex subset of M 2 . The dotted lines denote the 
smallest surrounding rectangle. 



a finite union of unit squares. Note that a polyomino is thus defined up to a translation as 
a subset of the plane IR 2 in which Z 2 is embedded. 

The area of a polyomino is the number of its squares. The boundary of a polyomino is the 
collection of unit edges of the dual lattice which belong only to one of the unit squares of the 
polyomino itself. The perimeter of a polyomino is the cardinality of its boundary, namely, 
the number of unit edges of the dual lattice which belong only to one of the unit squares 
of the polyomino itself. In other words the perimeter counts the number of interfaces on Z 2 
between the sites inside the polyomino and those outside. 

A polyomino is connected if and only if its interioj is a connected subset of M. 2 . A 
polyomino is convex if and only if its intersection with any line parallel to the coordinate 
axes of Z 2 is convex. 

A polyomino is monotone if and only if its perimeter is equal to that of its smallest 
surrounding rectangle. A polyomino is minimal if and only if any other polyomino with the 
same area has perimeter greater or equal to that of the polyomino itself. 

Lemma 6.16 The following holds: 

1. the perimeter of a connected polyomino is greater than or equal to that of its smallest 
surrounding rectangle; 

2. a connected polyomino is convex if and only if it is monotone; 

3. a minimal polyomino is connected and convex. 

The centers of the unit squares forming a connected polyomino are a nearest neighbor connected subset 
of Z 2 . As in [1] here a polyomino is not necessarily connected, different definitions, see for instance [6], can 
be found in the literature. 
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Proof of Lemma 16. 161 Call vertical and horizontal the two directions defined by the 
coordinate axes of Z 2 . Define, also, the notion of left, right, top, and bottom. Let a column 
(resp. row) be the subset of IR 2 obtained by considering the union of all the unit squares 
whose center belong to the same vertical (resp. horizontal) line. 

Item [1] since, by hypothesis, the polyomino is a connected subset of R 2 , each row and 
each column intersecting the smallest rectangle surrounding the polyomino intersects the 
boundary of the polyomino in at least two unit segments of the dual lattice. This remark 
yields the statement. 

Item [2j by definition of convex polyominoes it follows that a connected polyomino is 
convex if and only if each vertical and each horizontal line passing through the sites of Z 2 
intersects the polyomino in a convex not empty subset of the line itself. Thus, a polyomino 
is convex if and only if each vertical and each horizontal line passing through the sites of Z 2 
intersects the boundary of the polyomino in exactly two unit segments of the dual lattice. 
This remark yields the statement. 

Item |3j the fact that a minimal polyomino is connected is quite obvious. Indeed, it 
is sufficient to identify its maximal connected components and dispose them in a chain 
according to the following rule: excepted for the last component in the chain, the right part 
of the boundary of each component intersects the left part of the boundary of the component 
following it in the chain at least on a unit segment of the dual lattice. This construction 
yields a connected polyomino whose perimeter is smaller than the perimeter of the original 
polyomino by a quantity which at least equal to twice the number of maximal connected 
components minus one. 

We now prove that a minimal polyomino is convex. Consider a minimal polyomino c 
and, by absurdity, assume it is not convex. We shall construct another polyomino with the 
same area and smaller perimeter. This will yield the statement by contradiction. 

Since the polyomino is connected, by items Q] and [2] of this lemma we have that the 
perimeter P of the polyomino c is strictly greater than the perimeter R of its smallest 
surrounding rectangle r, that is to say P > R. 

As a main tool in the sequel of the proof we shall use horizontal and vertical projections as 
those introduced in [TJ Section 2, page 6], see also [191 Figure 8]. We first construct a second 
polyomino c\ by projecting c vertically to its bottom, that is we consider the polyomino c\ 
such that (i) the number of unit cell in each column is equal to that of the polyomino c, (ii) 
the intersection between c\ and any column is a convex polyomino, and (iii) the boundary of 
the polyomino obtained by intersecting c\ and any column intersects the bottom horizontal 
edge of r. We then construct a third polyomino C2 by projecting c\ horizontally to its left, 
that is we consider the polyomino C2 such that (i) the number of unit cell in each row is 
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equal to that of the polyomino ci, (ii) the intersection between C2 and any row is a convex 
polyomino, and (iii) the boundary of the polyomino obtained by intersecting c 2 and any row 
intersects the left vertical edge of r. 

We note, now, that c 2 is connected and convex and that the smallest rectangle r 2 sur- 
rounding C2 is a subset of r. Since P > R, the perimeter P 2 of c 2 is equal to the perimeter 
of r 2 , and r 2 C r, we have we have that P > R > P 2 . Which is an absurd since the area of 
c is equal to that of C2. □ 



Lemma 6.17 For any n positive integer there exists two positive integers s and k, with 
< k < s, such that either (i) n = s(s — 1) + k or (ii) n = s 2 + k. The set of polyominoes 
of area n and minimal perimeter contains a rectangle of side lengths s and s — 1 with a bar 
long k attached to one of its longest sides in the case (i) and a square of side length s with 
a bar long k attached to one of its sides in the case (ii). 

Proof of Lemma [6 . 1 7L The lemma is nothing but a simple restatement of Theorem 2.2 
in [1]. □ 



Corollary 6.18 For any n positive integer, the square of the perimeter of the polyominoes 
of area n is bounded from below by IQn. 

Proof of Corollary 16.181 In the case (i) of Lemma 16.171 we have that the perimeter of the 
polyomino is P = As. Since the polyomino is contained in a square of side length s, we have 
that n < s 2 . Hence P 2 = 16s 2 > IQn. 

In the case (ii) of Lemma \6. 171 we have that the perimeter of the polyomino is P = 4s + 2. 
Since the polyomino is contained in a rectangle of side lengths s and s + 1, we have that 
n < s 2 + s. Since 

P 2 = (4s + 2) 2 = 16(s 2 + s) + 4 > 16n + 4 
the proof of the corollary is completed. □ 



6.2. Metastable states of the Blume-Capel model 

In this subsection we prove the Lemmas 14.81 and 14.91 and the Theorems 14.101 and 14.111 which 
allow to identify the metastable states of the Blume-Capel model and to deduce the expo- 
nential estimate of the exit time. 

Proof of Lemma 14.81 We construct the following path. The path is defined by the 
following simple rule 
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Figure 6.5: Schematic representation of path introduced at the beginning of the proof of 
Lemma [4 .81 Only the first part of the path, the one going from d to few steps after V c , is 
shown. 



- start from d; 

- transform one minus spin in zero; 

- iteratively let the zero cluster grow until is reached by adding a zero protuberance 
to the longest side of a rectangle and then filling the slice with zeros; 

- transform one zero spin in plus; 

- iteratively let the plus cluster grow until u is reached by adding a plus protuberance 
to the longest side of a rectangle and then filling the slice with pluses. 

Note that the rectangles that are drawn during the above iterative procedure are squares or 
quasi-squares (rectangles whose side legths differ by one). 



The path, see figure I6.5[ starts at d and, passing through V c reaches 0, then passes 
through Q c and, finally, reaches u. This path can be divided into two parts U\ G f2(d, 0) 
connecting d to and Co>2 G f2(0, u) connecting to u. By direct inspection, it is very easy 
to show that 

$ W1 - H{d) = $ W2 - H(0) = r c 

and that U\ and U2 attain their maximal height respectively only in the states in V c and Q c . 
This yields the lemma. □ 

Now, we introduce the set C X as the set of states in X such that the number of 
minus spins is equal to 

|A| - [4(4 - 1) + 1] 



metadegenere.tex - 4 marzo 2013 



29 



21:47 



In the following lemma we characterize the set F{X d ) of the minima of the energy of X d . 



Lemma 6.19 Consider the zero chemical potential Blume-Capel model (fJ~2L) and assume 



that Condition \4- 7| is satisfied. Let a E X d and N a C A be the set of sites i <G A such that 
a(i) 7^—1. We have that 

1. the set N a is not a nearest neighbor connected subset of A winding around the torus A; 

2. if a G F(X d ) then a(i) = for any i e N a ; 

3. X d DV c ; 

I F{X d ) D V c ; 

5. H{F{X d )) = H{d) + T c . 



Proof of Lemma 16. 191 Item [T] recalling (14. 4p we have that 

|Ay = M4-l) + l<(f + l)f + l = | + f + l<i + ;| + i 

where in the last inequality we have used that, by Condition 14.7} we have h < 1. The item 
follows since, by Condition 14. 7\ we have that |A| > 49/ 7i 4 . 

Item[2j let a G X d and let r, £, m not negative integers such that r is the number of pluses 
in a, i is the number of plus-minus interfaces, and m is the number of plus-zero interfaces. 
Let a' be the configuration obtained by flipping to zero the pluses in a. By (14.21) . it follows 
that 

H(a') - H(a) = rh - 3£ - m < rh - (£ + m) 

Let C(a) be the polyomino obtained by collecting all the unit squares associated with 
a plus spin of a. Note that its perimeter is equal to £ + m and its area to r. Since the 
polyomino C(a) does not wind around the torus, we can apply Corollary 16.181 and obtain 
that {£ + m) 2 > 16r. This bound on the number of interfaces between the pluses and the 
other two types of spins implies that 

H(a')-H{a) < rh - Ay/r~ 

By studying the parabola of algebraic equation y = hx 2 — Ax it is immediate to show 
that H(a') — H(a) < provided y/r < A/h. This condition is easily proven to be true by 
using that r < £ c (£ c — 1) + 1. Indeed, by performing the same estimate as in the proof of 
item [TJ above, we have that 

r < 4(4 - 1) + K ^ 
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which yields the desired bound y/r < 4/h. 

We can finally conclude that H(a') < H(a), which completes the proof of the item. Note 
that the bound is not strict because, in principle, we could have no plus spin in a, that is to 
say r = 0. 

Item[3j since any configuration in V c has |A| — [£ c (£ c — 1) + 1] minus spins, it follows that 
V c C X d . 

Item HJ for any a G F(Xd) item |2] above states that the spins not equal to minus are 
equal to zero. Thus, F{X&) is a subset X^o C. X d made of those configurations such that all 
the spins that are not minus are zero. Consider a G X^ t o and denote by C(a) the polyomino 
obtained by collecting all the unit squares associated with a spin zero. Remarked that the 
area of the polyomino C(cr) is equal to [£ c (£ c — 1) + 1], by (14.21) we have that 

H{a) - H(d) = Perimeter (CV)) - h[£ c (£ c - 1) + 1] 

This quantity is minimal if and only if the perimeter of the polyomino C (a) is minimal. We 
can then conclude that the configurations in F(X&) are made of [£ c {£ c — 1) + 1] zeros in the 
sea of minuses and that the polyomino obtained by collecting all the unit squares associated 
with a zero spin has area [£ c {£ c — 1) + 1] and minimal perimeter. 

Since item [1] above ensures that for any a G X&q the polyomino C (a) does not wind 
around the torus, we can then apply Lemma [6.171 (case (i)) and conclude that V c C F{X^)- 

Item [5) immediate consequence of item H] above and the definition (14. 6 1) of T c . □ 

Proof of item[T]of Lemma [4.91 We have to compute <&(d, u) — H(d) and $(0, u) — H(0). 
In order to show that they are both equal to T c we show first that r c is an upper bound for 
both. This can be done easily, indeed by Lemma 14.81 it follows 

$(d, u) - H(d) < T c and $(0, u) - H(0) < T c (6.1) 

In order to get the lower bound the two problems have to be treated separately. We 
focus, first, on the problem $(d,u) and proceed as follows. 

Recall the definition of X^ given at the beginning of this subsection. Since the connec- 
tivity matrix q(a, rj) is different from zero only if the two configurations a and rj differ for 
the value of one single spin, we have that for any not negative integer n smaller than |A|, 
any path in Q(d, u) must necessarily visit the set of configurations with n minus spins. It 
then follows that any path in Q(d, u) must necessarily pass through X^. 

Item|5]in Lemma 16. 191 states that H(F(Xd)) = H(d) + T c . The above remarks imply that 
any path u G f2(d, u) is such that $ w > H(F(X d )) = H(d) + r c , which yields the desired 
lower bound $(d, u) - H(d) > T c . 
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The lower bound $(0,u) — H(0) > H(Q C ) — H(0) = T c can be achieved, "mutatis 
mutandis," with the same argument. We do not enter into the details, we just remark that 
the manifold must be replaced by Xq defined as the set of configurations in which all the 
spins are zeros excepted for £ c (£ c — 1) + 1 which are pluses. □ 

Before starting the proof of item [2] of Lemma 14.91 we need to state a technical lemma on 
the zero chemical potential Blume-Capel Hamiltonian ensuring that the energy (14. 2 p of a 
configuration is decreased if a minus spin having at most two minuses between its nearest 
neighbors or a minus spin having three neighboring minus and a neighboring plus is flipped 
to zero. Moreover, the lemma states that the energy of a configuration is decreased if a zero 
spin with at least two pluses and no minus among its nearest neighbors is flipped to plus. 

For any a G X, i G A, and a G { — 1,0, +1} we let a l,a be the configuration such that 
crha (j) — f° r a ll J £ A and j ^ i and a l ' a (i) = a. Note that a %,a differs from a by at 
most the value of the spin associated with the site %. 



Lemma 6.20 Consider the zero chemical potential Blume-Capel model ( 4-2 ) and assume 
that Condition ^. 7| is satisfied. Let a G X 

1. if there exists a site i such that a(i) = — 1, three of the nearest neighbors of i have 
associated spin equal to —1, and the fourth has associated spin equal to +1, then 



H(a)-H(a 



i,0\ 



2. if there exists a site i such that a(i) = — 1 and there are at most two nearest neighbor 
of % such that the associated spin in a is —1, then 

H(a) - H(a l '°) > h 

3. if there exists a site i such that a(i) = 0, there are at most two nearest neighbors 
of i such that the associated spins are 0, and the remaining nearest neighbors have 
associated spin equal to +1, then 

H(a) - H(a l ' +1 ) > h 



Proof of Lemma [6.201 Items [T] and [2) let n+,n_,no > be the number of plus, minus, 
and zero spins, respectively, among the four nearest neighbors of the site i. Note that 
n + + n_ + no = 4. From (14.21) . use also that here a(i) = —1, it follows immediately that 

H(a) - H(a ifl ) = 4n + + n + h - [n_ + n+] = h + 3n + - n_ + n 
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By using that n + + no = 4 — n_ we get 

H{a) - H (a ifi ) = h + A + 2(n + - n_) 

Item [1] follows by the formula above with = 3 and n + — 1; item [2] follows by the formula 
above by noticing that n_ < 2 and n + > 0. 

Item[3j let n+, no > be the number of plus and zero spins, respectively, among the four 
nearest neighbors of the site i. Note that, by hypothesis, n + + n = 4. From (14.21) . use also 
that here a(i) = 0, it follows immediately that 

H(a) - H(a l > +1 ) = n + - [n - h] = 4 - 2n + h 

where we have used n + = 4 — uq. The item finally follows by noticing that, by hypothesis, 
n < 2. □ 



Proof of item [2] of Lemma 14.91 Let a G X \ {d, 0, u}, we have to prove that V a < T c ; 
that is we have to find a configuration at energy smaller than H(a) and a path connecting 
a to such a configuration such that its height is smaller than H(a) + T c . 

We first assume that a is such that at least one spin is equal to minus one and we 
distinguish several different sub-cases. 

Case 1. There exists a plus-minus interface in a, that is there exist i,j G A nearest neighbors 
such that a(i) = — 1 and a(j) = +1. From items [1] and [2] in Lemma 16.201 we have that 
H(a) - H(a l '°) > h. It then follows H(a) > H{a ifi ) and, hence, V a = 0. 

Case 2. No plus-minus interface exist in a and there exist a site % G A such that a(i) = —1 and 
at least two of its nearest neighbors are occupied by spin zero. From item [2] in Lemma 16.201 
we have that H(a) - H(a ifi ) > h. It then follows H{a) > H(a ifi ) and, hence, V a = 0. 

Case 3. No plus-minus interface exist in a and for any i e A such that a(i) = —1 there 
is at most one nearest neighbor occupied by spin zero. Let 7 be the collection of the unit 
segments of the dual of A separating two sites with which are associated a minus and a 
zero spin; call those sites adjacent to 7. This collection 7 is made of maximal connected 
components which are either rectangles or straight annuli winded around the torus. Note 
that by the characterization of a given in this case 3 we have that each minus associated 
with a site adjacent to 7 has one single zero (the one associated to the neighbor adjacent to 
7) among its nearest neighbors. 

Case 3.1. Suppose that one of the connected component of 7 is a straight annulus a winding 
around the torus. The annulus separates a stripe of zeros by a stripe of minuses (both 
winding around the torus). Since each minus associated with a site adjacent to a has one 
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single zero among its nearest neighbors, it follows that the minus stripe occupies at least 
two lattice rows or columns. Then, let a' be the configuration obtained by flipping to zero 
all the minuses of a associated with sites adjacent to a and note that, by ( 14. 2ft . 

H(a') - H{a) = -hL H(a') < H(a) 

where we recall that L is the side length of A. Moreover, consider the L + 1 long path 
connecting a to a' constructed by flipping to zero first one of the minus spin associated with 
a site adjacent to a and then the remaining L — 1 minuses so that at each step one minus 
with two zeros among its nearest neighbors is flipped. The height of this path is H (a) + 2 — h. 
Since r c > 2 — h we have that V a <T C . 

Case 3.2. Suppose that none of the connected component of 7 is a straight annulus winding 
around the torus and assume that there exists a rectangular component such that one of its 
four sides, say p, has length \p\ larger or equal to £ c . Let a' be the configuration obtained 
by flipping to zero all the minuses of a associated with sites adjacent to p. Since each minus 
associated with a site adjacent to p has one single zero among its nearest neighbors, by (14. 2p . 
we have that 

H(a') - H(a) = 2- h\p\ < 2 - hl c < 2 - h- =>• H(a') < H(a) 

h 

where we have used that £ c > 2/h (recall (14. 3p ). Moreover, consider the \p\ + 1 long path 
connecting a to a' constructed by flipping to zero first one of the minus spin associated with 
a site adjacent to p and then the remaining \p\ — 1 minuses so that at each step one minus 
with two zeros among its nearest neighbors is flipped. The height of this path is H(a) + 2 — h. 
Since r c > 2 — h we have that V a < T c . 

Case 3.3. Suppose that none of the connected component of 7 is a straight annulus winding 
around the torus and assume that there exists a rectangular component such that one of 
its four sides, say p, has length \p\ = 1. Let a' be the configuration obtained by flipping 
to minus the zero of a associated with sites adjacent to p. By considering the two possible 
situations, that is to say the lxl square and the 1 x i rectangle, by (" 14 . 2 [) we have that 

H(a') - H(a) < max{-4 + h, -2 + h} = -2 + h< => H(a') < H(a) 

where we have used that h < 1. Since a is transformed in a' in one step, we have that 

v a = < r c . 

Case 3.4. Suppose that none of the connected component of 7 is a straight annulus winding 
around the torus and assume that there exists a rectangular component such that one of 
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its four sides, say p, has length \p\ = 2. Let a' be the configuration obtained by flipping to 
minus the two zeros of a associated with sites adjacent to p. By (14. 2 p we have that 

H(a') - H{a) = -2 + 2h < H{a') < H{a) 

where we have used that h < 1. Moreover, consider the path long three connecting a to a' 
constructed by flipping to minus the two zero spins associated with the sites adjacent to p. 
The height of this path is H(a) + h. Since T c > h we have that V a < T c . 

Case 3.5. Suppose that none of the connected component of 7 is a straight annulus winding 
around the torus and assume that all the rectangular components have side lengths larger 
or equal to three and smaller or equal to £ c — 1. Note that this case can be empty if £ c = 3. 
Consider one of this rectangular component and let p be one of its four sides; note that 
3 < \p\ < £ c — 1- For the sake of clearness assume that p is vertical. 

Case 3.5.1. Suppose that the column of sites at distance one from those occupied by zeros and 
adjacent to p (in other less precise words the "second" column inside the rectangle starting 
from p) is occupied only by spin zero. Let a' be the configuration obtained by flipping to 
minus the |p| zeros of a associated with sites adjacent to p. By (14. 2 h we have that 

H(a') - H(a) = -2 + h\p\ < -2 + h{£ c - 1) < -2 + h- < => H(a') < H(a) 

where we have used that \p\ < £ c — 1. Moreover, consider the \p\ + 1 long path connecting a 
to a' constructed by flipping to minus, one after the others, the zero spins associated with 
the sites adjacent to p and by flipping, in the first \p\ — 1 steps, one of the two zeros with 
two minuses among its four nearest neighbors; in the final step the remaining spin zero is 
finally flipped to minus. The height of this path is H(a) + h(\p\ — 1). Since T c > h(£ c — 2) 
we have that V a < L c . 

Case 3.5.2. Suppose that the column of sites at distance one from those occupied by zeros and 
adjacent to p (in other less precise words the "second" column inside the rectangle starting 
from p) is occupied at least by one spin plus. Let a' be the configuration obtained by flipping 
to zero the pluses of a associated with sites at distance one from those adjacent to p and 
occupied by zeros. 

Case 3.5.2.1. Suppose \p\ — 2 = 1. The configuration a is transformed into a' by flipping to 
zero a plus spin having at least three zeros among its nearest neighbors and no minus. It 
then follows that 

H(a') - H(a) < max{-4 + h, -2 + h} = -2 + h < H(a') < H(a) 
where we have used that h < 1. Since a is transformed into o' in one step we have that 

K < < r c . 
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Case 3.5.2.2. Suppose \p\ — 2 > 2. The configuration a is transformed into a' by flipping to 
zero at most \p\ — 3 plus spin having two zeros and two pluses as nearest neighbors and at 
least one plus spin having no minus and at most one plus among its nearest neighbors. It 
then follows that 

H(a') - H(a) < max{-4 + h, -2 + h} + (\p\ - 3)h < -2 + (|p| - 2)h 
Hence, recalling that \p\ < £ c — 1, we have 

H(a') - H(a) < -2 + (4 - 3)h < -2 + | h - 2h < -2h H(a') < H(a) 

Moreover, consider the path connecting a to a' constructed by flipping to zero, one after 
the others, the plus spins with two neighboring pluses and then those with at most one 
neighboring plus. The height of this path is smaller that H(a)+h(\p\ — 3). Since T c > h(l c —A) 
we have that V a < T c . 

In order to complete the proof, we finally have to consider those configurations a G 
X\ {d, 0, u} such that none of the spins is minus. The proof can be achieved with arguments 
very similar to those developed above; we just sketch the idea. First of all we consider the 
collection of the unit segments of the dual of A separating two sites with which are associated 
a zero and a plus spin. From item [3] in Lemma 16.201 it follows that this collection is made of 
maximal connected components which are either rectangles or straight annuli winded around 
the torus. We then have to consider all the cases analogous to those taken into account above. 

□ 

Proof of Theorem 14.101 The theorem follows from Lemma 14.91 and Theorem 12.41 □ 

Proof of Theorem 14.111 The theorem follows from Theorem 14.101 above and [22J, Theo- 
rems 4.1 and 4.9]. □ 

6.3. Escape mechanism 

In this subsection we prove Lemma T4. 121 and Theorem 14.131 allowing the identification of the 
escape mechanism as the nucleation of the critical droplet. 

Proof of Lemma f4. 121 Item[T] the proof of this item is divided into two parts. In the first 
part we characterize the set F{X^) of the minima of the energy of the set C X defined 
at the beginning of Subsection 16.21 In the second part we show that any path connecting d 
to u must necessarily pass through V c . 
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Figure 6.6: Typical configuration of Fi(X d ) (left) and F 2 (A"d) (center and right). Zeros are 
associated with the sites inside the solid lines and minuses outside. 

First part of the proof. From items Q] and [2] of Lemma 16.191 it follows that for any a G 
F(X d ) the polyomino C(cr) associated with the zero spins does not wind around the torus. We 
can then apply item [3] in Lemma [6.161 and deduce that C(a) is convex. Thus, we partition 
the set F{X<±) into two disjoint subsets F-\_{X&) and F 2 (X < i) defined as follows (see, also, 
figure WM : the set Fi(X d ) is the set of configurations a G F(X d ) such that the boundary of 
the polyomino C (a) intersects each side of the boundary of its smallest surrounding rectangle 
on a set of the dual lattice Z 2 + (1/2, 1/2) made by at least two pairwise consecutive unit 
segments. The set ^(^d) is the set of configurations a G F(X d ) such that the boundary of 
the polyomino C(a) intersects at least one of the four sides of the boundary of its smallest 
surrounding rectangle on a set made of a single unit segments. Note that, since C(a) is 
convex, there cannot exist multiple unit protuberances intersecting the same side of the 
smallest surrounding rectangle. 

In figure IBTBl the two configurations in i^O^d) have been represented by a rectangle of 
zeros with a unit zero protuberance placed either on the longest (picture on the right) or on 
the shortest side (center picture). In principle other situations should be taken into account, 
indeed the boundary of the polyomino associated with the zero component of the considered 
configuration could intersect the other three sides of the boundary of its smallest surrounding 
rectangle in proper subsets of the side itself. In other words the rectangle on which the unit 
protuberance is placed could be not fully occupied by zeros. 

We prove, now, that the two situations depicted in the picture indeed cover all the 
possible cases. Consider a G Fi{X&). Let l c + k and l c + s be the side lengths of the smallest 
rectangle surrounding the polyomino C(er) associated with a with k, s G Z. Assume, also, 
that a single protuberance U is placed on the side of length £ c + s. Let R(cr) be the smallest 
rectangle surrounding the polyomino obtained by removing the unit cell U from C(a); note 
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that R(cr) is a rectangle with side lengths 4 + s and £ c + k — 1. Since the polyomino C(er) 
is convex, we have that its perimeter is equal to that of the smallest surrounding rectangle, 
namely, 

Perimeter(C(a)) = 2(4 + k) + 2(4 + s) = 44 + 2(k + s) 

Moreover, from item H] in Lemma 16.191 it follows that Perimeter (C (a)) = 44 and, hence, 
k + s = 0. 

Note, now, that the area of the rectangle R(cr) is 

(4 + s)(4 + k - l) = 4(4 - l) + (Jfe + s)4 + &s - s = 4(4 - l) - s 2 - s 

where we have used k + s = 0. Since the area of the polyomino C(er) is 4(4 ~ 1) + 1, the 
number of zeros inside R(cr) is 4(4 — 1)- Hence the area of R(cr) must be larger or equal to 
4(4 — 1); it then follows that 

4(4 - i) - s 2 - s > 4(4 - l) -s 2 - s > o 

Since s is an integer, the above inequality is satisfied only for s = and s = — 1; note that 
in both cases the inequality is indeed an equality We can then give a full characterization 
of the set -F 2 (Ad). Indeed we can write 

with i^aC^d) = 'Pc and i^&^d) the set of configurations of F(X&) such that the zeros 
occupy a rectangle of side lengths 4 — 1 and 4 and a unit protuberance placed on one of 
the two shortest sides. 

Second part of the proof. We prove, now, that any path connecting d to u necessarily 
pass through V c . The proof is organized in several steps. 

Step 0: for any path oj G fl(d, u) there exists a positive integer i such that uji G X d . See 
the proof (lower bound) of item [1] of Lemma 14.91 

Step 1: for any path u G f2(d, u) such that $ w — H(d) = T c there exists a positive integer 
i such that Ui G F(X&). This property follows from step and item H] in Lemma [6. 191 

Step 2: for any path uj G f2(d, u) such that $ w — H(d) = T c we denote by f(oj) the set 
of positive integers such that for any i G f{oo) one has that Ui G F(X&) and is made of 
I A | — 4(4 — 1) minus spins and 4(4 — 1) zeros. From step 1 above it follows that f(oS) is 
not empty. 

Step 3: for any path u G f2(d, u) such that ^ — H(d) = T c one has that Ui G -^(^d) 
for any % G f(oo). Indeed, assume by absurdity that that there exists % G f{w) such that 
Ui G F 1 (X d ). By definition of f{ui) the configuration Ui is transformed into by flipping 
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Figure 6.7: Paths constructed in the second part of the proof of Lemma [4.121 

to minus one of the zeros in Ui. From the definition of Fi(Xd) it follows that all the zero 
spins in Ui have at least two zeros among their nearest neighbors. Then by (14. 2 j) it follows 
easily that H(u>i-i) > H{uJi). Hence, 

^ > H{ui- X ) > H(ui) = H(F(X d )) = H(d) + T c 

which is an absurd. We can then conclude that Ui G -^('fd) for any % G f(ui). 

Step 4: for any path oj G Q(d, u) such that — H(d) = T c , if i is a positive integer such 
that Ui G F(X d ) then Ui_i,u i+1 £ X d . Indeed, assume by absurdity that G X d . Then 
Ui+i is necessarily obtained by Ui by flipping to plus a zero spin, otherwise the number of 
minus spins would be changed (recall that all the configurations in X d have the same number 
of minuses). Since any zero spin of has no plus among its nearest neighbors, from (14. 2 \ it 
follows that H(u i+ i) > H(ujj). Hence, 

$ w > H{oo i+l ) > H{ Ui ) = H{F{X A )) = H{d) + T c 

which is an absurd. We can then conclude that uji + i G" X d . Similarly we also prove that 
0Ji-i G" X a . 

Step 5: consider a path u G Q(d, u) such that $ w — H(d) = T c and assume, by absurdity, 
that it does not pass through V c , that is to say uj% G^ V c for any integer i. 

From step 4 (see figure 16.71) it follows that there exists n not consecutive integers i\ < 
%2 < ■ ■ • < i n -i < i n such that u ik G F(X d ), for k = 1, . . . , n. If n > 2, each sub-path 
(cu ik+ i, . . . , Wj fe+1 _i), for k = 1, . . . , n — 1, belongs either to the subset X d < of X made of all 
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the configurations with number of minus spins smaller or equal to |A| — [4(4 — 1) + 1] — 1 
or to the subset X^> of X made of all the configurations with number of minus spins larger 
or equal to |A| - [4(4 - 1) + 1] + 1. 

Let j be the smallest integer in {1, . . . , n} such that the path (o;^, O0i +i, . . . , u) does not 
visit X dt> - Since i G X^y, from step 3 and from the absurd hypothesis it follows that 

^ e F 2>b (X d ). 

Since H^u^) = H(d) + r c , starting from the path must necessarily decrease the 
energy. Since it is not possible to flip to minus the unit zero protuberance (otherwise the 
path would enter Xd,>), the only move that decreases the energy is flipping to zero a minus 
with two zeros among its nearest neighbors. So that H(u)i j+ i) = H(d) + T C — h. 

Starting from oj^+i only moves which increase the energy of at most h are allowed. So 
that the sole possible moves are: flipping to zero a minus with two zeros among its nearest 
neighbors or flipping to minus a zero with two minuses among its nearest neighbors. 

With these types of moves only configurations a such that the polyomino C(a) associated 
with the zero spin component is convex and its smallest surrounding rectangle has side 
lengths 4 + 1 and £ c — l. In this set the smallest energy is that of the configuration in which 
C(cr) is precisely the rectangle with side lengths 4 + 1 and £ c — 1- 

Besides the already mentioned moves, the one to which competes the smallest energy 
increase 2 — h is flipping to zero a minus with one zero and three minuses among its nearest 
neighbors. Thus, we have that 

& u ~ H(d) > 2(4 + 1) + 2(4 - 1) - h(£ c + 1)(4 - 1) + 2 - h 

Hence 

- H(d) > 44 - M4(4 - 1) + 1] - h(£ c - 1) + 2 > r c 

where, in the last inequality, we have used (14. 5ft . (14. 6ft . and the bound 

2 

2 - h(£ c - 1) > 2 - h- > 

where we recalled the upper bound (14.41) . 

We finally got an absurd. We then have that the path uj has to visit V c . This completes 
the proof of item [1] of the lemma. 

Item [2) the proof can be achieved, "mutatis mutandis," with the same arguments used 
in the proof of item [TJ We do not enter into the details, we just remark that the manifold 
Xd must be replaced by X defined as the set of configurations in which all the spins are 
zeros excepted for 4(4 — 1) + 1 which are pluses. □ 

Proof of Theorem 14.131 In the proof of this theorem we use results in [22] , in particular 
the definition of gate on page 603 and the Theorem 5.4. The Lemma T4. 121 in Section H] above 
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implies that V c is a gate for the pair of configurations d and u and Q c is a gate for the pair 
of configurations and u. Then the theorem follows from Theorem 5.4 in [22]. □ 

Proof of Theorem 14. 141 As stated in item [1] of Lemma 14.81 (see also figure I6.5p . there 
exists a path joining d to u passing through any state in V c and attaining its maximal height 
only in this point. This implies that any subset of X not containing V c is not a gate, in 
particular it is not minimal. 

The second statement of the theorem can be prove similarly. □ 

References 

[1] L. Alonso, R. Cerf, "The three dimensional polyominoes of minimal area." The Elec- 
tronic Journal of Combinatorics 3, #R27 (1996). 

[2] J. Beltran, C. Landim, "Tunneling and Metastability of Continuous Time Markov 
Chains." Journ. Stat. Phys. 140, 1065-1114 (2010). 

[3] S. Bigelis, E.N.M. Cirillo, J.L. Lebowitz, E.R. Speer, "Critical droplets in metastable 
probabilistic cellular automata," Phys. Rev. E 59, 3935 (1999). 

[4] M. Blume, "Theory of the First-Order Magnetic Phase Change in UO2." Phys. Rev. 
141, 517 (1966). 

[5] M. Blume, V.J. Emery, R.B. Griffiths, "Ising Model for the A Transition and Phase 
Separation in He 3 -He 4 Mixtures." Phys. Rev. A 4, 1071-1077 (1971). 

[6] M. Bousquet-Melou, J.M. Fedou, "The generating function of convex polyominoes: the 
resolution of a g-differential system." Discrete Mathematics 137, 53-75 (1995). 

[7] A. Bovier, "Metastability: a potential theoretic approach." Proceedings of ICM 2006, 
EMS Publishing House, 2006, pp. 499-518. 

[8] A. Bovier, M. Eckhoff, V. Gayrard, M. Klein, "Metastability and low lying spectra in 
reversible Markov chains." Comm. Math. Phys. 228, 219-255 (2002). 

[9] A. Bovier, F. den Hollander, and F.R. Nardi, "Sharp asymptotics for Kawasaki dynamics 
on a finite box with open boundary." Probab. Theory Relat. Fields 135, 265-310 (2006). 

[10] A. Bovier, F. Manzo, "Metastability in Glauber dynamics in the low-temperature limit: 
beyond exponential asymptotics." Journ. Stat. Phys. 107, 757-779 (2002). 



metadegenere.tex - 4 marzo 2013 



11 



21:47 



[11] H.W. Capel, "On possibility of first-order phase transitions in ising systems of triplet 
ions with zero-field splitting." Physica 32, 966 (1966). 

[12] E.N.M. Cirillo, F.R. Nardi, "Metastability for the Ising model with a parallel dynamics." 
Journ. Stat. Phys. 110, 183-217 (2003). 

[13] E.N.M. Cirillo, F.R. Nardi, A.D. Polosa, "Magnetic order in the Ising model with parallel 
dynamics." Phys. Rev. E 64, 57103 (2001). 

[14] E.N.M. Cirillo, F.R. Nardi, C. Spitoni, "Metastability for reversible probabilistic cellular 
automata with self-interaction." Journ. Stat. Phys. 132, 431-471 (2008). 

[15] E.N.M. Cirillo, F.R. Nardi, C. Spitoni, "Competitive nucleation in reversible Proba- 
bilistic Cellular Automata." Phys. Rev. E 78, 040601 (2008). 

[16] E.N.M. Cirillo, F.R. Nardi, C. Spitoni, "Competitive nucleation in metastable systems." 
Applied and Industrial Mathematics in Italy III, Series on Advances in Mathematics 
for Applied Sciences, Vol. 82, Pages 208-219, 2010. 

[17] E.N.M. Cirillo, E. Olivieri, "Metastability and nucleation for the Blume-Capel model. 
Different mechanisms of transition," Journ. Stat. Phys. 83, 473-554 (1996). 

[18] R.J. Glauber, "Time-Dependent Statistics of the Ising Model." J. Math. Phys. 4, 294 
(1963). 

[19] F. den Hollander, F.R. Nardi, A. Troiani, "Kawasaki Dynamics with Two Types of 
Particles: Stable/Metastable Configurations and Communication Heights." Journ. Stat. 
Phys. 145, 1423-1457 (2011). 

[20] F. den Hollander, F.R. Nardi, A. Troiani, " Metastability for Kawasaki dynamics at low 
temperature with two types of particles." Electronic Journal of Probability 17, 1-26 
(2012). 

[21] G. Grinstein, C. Jayaprakash, Y. He, "Statistical Mechanics of Probabilistic Cellular 
Automata." Phys. Rev. Lett. 55, 2527-2530 (1985). 

[22] F. Manzo, F.R. Nardi, E. Olivieri, E. Scoppola, "On the essential features of metasta- 
bility: tunnelling time and critical configurations." Journ. Stat. Phys. 115, 591-642 
(2004). 

[23] F. Manzo, E. Olivieri, "Dynamical Blume-Capel model: Competing metastable states 
at infinite volume." Journ. Stat. Phys. 104, 1029-1090 (2001). 



metadegenere.tex - 4 marzo 2013 



12 



21:47 



[24] F.R. Nardi, C. Spitoni, "Sharp Asymptotics for Stochastic Dynamics with Parallel Up- 
dating Rule." Journ. Stat. Phys. 146, 701-718 (2012). 

[25] E. Olivieri, E. Scoppola, "Markov chains with exponentially small transition probabil- 
ities: First exit problem from a general domain. I. The reversible case," Journ. Stat. 
Phys. 79, 613-647 (1995). 

[26] D.J. Wales, Energy landscapes. Cambridge University Press, Cambridge, UK, 2003. 



metadegenere.tex - 4 rnarzo 2013 



13 



21:47 



